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Inelastic x-ray scattering (IXS) now a widely used technique for studying the dynamics of electrons 
in condensed matter. We previously posed a solution to the phase problem for IXS [P. Abbamonte, 
et. al., Phys. Rev. Lett. 92, 237401 (2004)], that allows explicit reconstruction of the density 
propagator of a system. The propagator represents, physically, the response of the system to an 
idealized, point perturbation, so provides direct, real-time images of electron motion with attosecond 
time resolution and A-scale spatial resolution. Here we show that the images generated by our 
procedure, as it was originally posed, are spatial averages over all source locations. Within an 
idealized, atomic-like model, we show that in most cases a simple relationship to the complete, un- 
averaged response can still be determined. We illustrate this concept for recent IXS measurements 
of single crystal graphite. 

PACS numbers: ??, ?? 



CO 



I 

C 

O 

o 



> 

in 

On 
l> 
O 

o 

OS 

o 



I. INTRODUCTION 

Attoscience is an emerging discipline that exploits re- 
markable achievements in HHG laser technology to study 
the dynamics of electrons in atoms, molecules, and con- 
densed matter with sub- femtosecond time resolution. [l|, 
ULEI These advances have been paralleled by also remark- 
able achievements in technology for x-ray experiments, 
both at the source and the endstation level. In partic- 
ular, new approaches to inelastic x-ray scattering (IXS) 
provide a kinematic flexibility not previously accesible to 
energy-loss scattering techniques. [4[ 

We previously If the scattering cross section is sampled 
over a sufficiently broad range of momentum and energy, 
it can be inverted to explicitly reconstruct the density- 
density propagator of a material. The propagator, phys- 
ically, represents the density induced in a medium by an 
idealized, point perturbation. Inversion therefore pro- 
vides direct, real time images of electron dynamics, with 
resolutions exceeding 20 as in time and 1 A in space. [5[ 
We have used this method, for example, to show that the 
valence exciton in LiF is a Frenkel exciton.Q 

This inversion method, as it was originally posed @, 
has a limitation: it produces spatially averaged images. 
To understand the origin of this limitation, recall that 
the (retarded) density propagator of a system is defined 



x (x 1 ,x 2 ,t) = -i/h<g\ [h(xx,t),n(x 2 ,0)] \g > 6{t) (I) 

where \g > represents the ground state, and h(x,t) is the 
density operator 



h{x,t) = $ (x, t)ip[x, t). 



(2) 



Eq. I represents the amplitude that a point disturbance 
in the electron density at location x 2 at time t = will 
propagate to point xi at some later time t. 



The two spatial coordinates, x% and x%, are distinct. 
That is, the form of the response depends, in general, 
on the precise location of the source, (on top of an 
atom versus between atoms, for example). The Fourier 
transform of the propagator is therefore a function of two 
momenta 



X = X( k i,k 2 ,uj). 



(3) 



Unfortunately, in scattering there is only one momen- 
tum transfer, k. IXS experiments do not measure all 
parts of the response in eq. 3, but only its longitudinal 
(or "diagonal") components [7[ 



x(k,uj) = x{k,-k,w). 



(4) 



This restricts the amount of information available with 
the IXS imaging approach. 

If the system is translationally invariant, this restric- 
tion poses no limits on the method. In a homogeneous 
system the propagator depends only on the difference 
x = x\ — X2, i.e. 



X(x 1 ,x 2 ,t) = x(xi - x 2 ,t) = x{x,t), (5) 
its Fourier transform is a function of only one momentum 



X = X(k,w) 



(0) 



and IXS can provide a complete parameterization of the 
density propagator. 

All real systems are made of atoms, however, so are 
never truly homogeneous. One can, nonetheless, carry 
out the inversion procedure described in refs. 0, Q and 
construct detailed images of electron motion. But these 
images cannot be complete. Because they are based on a 
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FIG. 1: Periodic array of potential wells. For suitable choice 
of well depth and spacing, this system can be described in 
the tight binding approximation as a set of coupled harmonic 
oscillators. 
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FIG. 2: Contour plot of the prefactor 4>o(xi)il>o{x2) ■ 



subset of all Fourier components of the response, the im- 
ages must represent some kind of average of the complete 
response. For IXS inversion to be useful, it is necessary 
to articulate exactly what kind of average this is. 

In this paper we show that the images produced in 
an IXS inversion correspond to the complete response, 
x(x%,X2,t), averaged over all source locations, x 2 - In a 
homogeneous system, in which the response is indepen- 
dent of x 2 , the images can be considered complete. Even 
in an inhomogeneous system, however, a simple physical 
picture is still applicable. We show, within a simple but 
quite general model, that the overall size of the response 
scales with the electron density at the source location, 
x 2 . The average will, therefore, be dominated by source 
points of high electron density. We will show that, in 
cases of physical interest, the images will resemble closely 
the full response, x(. x ii x 2,i), hi which x 2 coincides with 
the peak in the electron density. We discuss this conclu- 
sion in the context of some recent IXS studies of single 
crystal graphite. 

II. INVERSION OF IXS DATA 

We begin by reviewing the mechanics of an IXS in- 
version. Inelastic x-ray scattering measures the dynamic 
structure factor S{k,oj), defined as 

S(k,u) = ^2\< f\n(k)\i >| 2 PiSiu-ujf + Ui) (7) 

where Pi — e~@ u)i /Z is the Boltzmann factor (/3 = h/kT). 
The Fourier transform of S is a space-time correlation 
function for the density. This correlation function, while 
interesting, is not causal so contains no direct information 
about dynamics. 

S(k,u>) is, however, related to the charge propagator 
through the fluctuation-dissipation theorem 



S{k,u) = - l _ e _ (j J m[ X {k,u)}. (8) 

The propagator is a truly causal quantity, so provides 
direct information about electron dynamics. 

Inverting this cross section, to reconstruct the propa- 
gator, requires four steps. The first is symmetrizing to 
convert from S to \. This is best done via the relation- 
ship @ 

Im[x[k, uj)} = -it [S(k, u) - S(k, -u)] (9) 

which does not require explicit specification of the tem- 
perature. Next, one must analytically continue the data. 
This can be done with simple, linear interpolationjg, 
Third, one performs a sine transform of this imaginary 
part 

X (k,t) = J^Im[x(kM]s^t (10) 

where it is understood that x(k, t) is zero for t < 0. Eq. 
(10) follows directly from the Kramers Kronig relations 
and is the step in the procedure in which the phase is de- 
termined and causality is imposed. Finally, one performs 
a normal, Fourier transform of the spatial coordinates to 
get x{x,t). 

This procedure has been used to image electronic mo- 
tion in liquid water to show that the exciton in LiF is 
a Frenkel exciton Q, and to study screening processes in 
graphite and grapheneQ. 

III. IMPLICIT SPATIAL AVERAGING 

We now stop to ask the following question: How do we 
interpret the output of this procedure for a system that 
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7(2:1,2:2) = 8(xi - x 2 ). 



(14) 



x(x) is the Fourier transform of the product 
g(ki, k 2 )x(ki, k 2 ), which means it can be written 
as the convolution integral 



X 2 ) 



) = / dx[dx 2 g(x' x ,x' 2 ) xipi - x[,x 2 - x' 2 ) 



(15) 



(16) 



Redefining the origin and writing in terms of x = x\ — x 2 
we arrive at 



FIG. 3: Plots of the function f n (t) for values of n out 
to fourth neighbors. Energy rapidly spreads from the first 
Wigner-Seitz cell (n — 0) into neighboring cells. 

does not have translational symmetry? If we carry out 
the procedure described in Section 2, we will be left with 
a function of a single spatial variable, x( x i t), which looks 
like the response of a homogeneous system. How should 
this function be interpreted? 

To answer this question, we must determine an ex- 
plicit relationship between this \(x,t) and the complete 
response function x ( x i> x 2i t). First, we write down the 
Fourier transform 

X(xt,x 2 ) = J ^^ x{kl ,k 2 )e ik ^+ ik ^, (11) 

where the time dependence has been suppressed. To sim- 
plify the discussion we have assumed that our system is 
one-dimensional, i.e. the coordinates x, x±, and 2:2 are 
scalars. 

Next, we recognize that if we insert the expression 
g(ki, k 2 ) — 27T(5(fci + k 2 ) into the integrand, 

X(xi,x 2 )= J ^l^lg(k 1 ,k 2 )x(k ll k 2 )e^ + ^ 

(12) 

cq. 11 reduces to 

X(xi - 2:2) = / ^ X(h,-h) e '*x(xx-«,) j (13) 

which we readily recognize to be the same as eqs. 4-5. 
This means that one can think of xi x i t) as the full func- 
tion x( x i , 2:2) having been Fourier filtered by the function 
g(ki,k 2 ). 

A more physical interpretation can be achieved by ap- 
plying the convolution theorem. The Fourier transform 
of g is given by 



X(x) = I dx' 2 x(x + x' 2l x' 2 ). (17) 

So, we see that the physical meaning of x{ x i t) is quite 
clear. It is simply the full susceptibility, averaged over 
all possible source locations, 2:2- No such average is done 
in time, however, so xi x , t) is still a direct probe of elec- 
tron dynamics. Further, even in an inhomogeneous sys- 
tem, xi x ,t) provides direct insight into the full response 
x(xi, x 2 , t). To illustrate why, it is useful to consider a 
specific case. 

IV. ONE ELECTRON IN A PERIODIC 
POTENTIAL 

To illustrate the value of x(x,t) for an inhomogeneous 
system, we consider the idealized case of a single electron 
traveling in a one-dimensional periodic array of harmonic 
wells, 

V{x) = ^ j v{x-na) (18) 

n 

where each well is described by 

v(x) = k{x 2 - r 2 )9(r - \x\), (19) 

a being the distance between wells and r the radius of 
each well, as shown in Fig. 1. For sufficiently large val- 
ues of k, the ground state of an isolated well will be a 
Gaussian 

(j)(x) = V2naexp~x 2 /2a 2 (20) 

with energy wo- 
lf the wells are weakly coupled the system can be de- 
scribed by a tight binding model with hopping parameter 



4 



w 

T3 
C 

o 

o 
o 
cn 
o 



») 

T3 

c 
o 

o 
o 
w 
o 



5000 
4000 




(a) 


3000 






2000 






1000 



• • 





5000 
4000 
3000 
2000 
1000 


0.25 
0.20 
0.15 
0.10 
0.05 
0.00 
-0.05 



am 





— x 2 = a/2 (magnified 10x) t \ 





! -1 1 2 

x/lattice parameters) 

FIG. 4: (a) Contour plot of x( x i, x 2,t) for X2 = 0, i.e. for the 
source located in the center of an "atom". The disturbance 
created is symmetric around the source point and propagates 
at the average group velocity of the band, (b) Same plot, but 
with X2 = a/2, i.e. for the source half way in between two 
atoms. The disturbance is now asymmetric, and no longer 
oscillates with a well-defined period, (c) Sections through the 
plots in (a) and (b) at t = lOOOas on a linear plot, showing 
that the magnitude of the response in (a) is larger by ~ 100 
times. 



th (to distinguish it from the time t). The states will be- 
come Bloch waves 



iP k {x) = Y,<t>{x-na)e ikna , 



(21) 



and the energy level will spread into a band with disper- 
sion 



oj(k) = 2t/j[l — cos(fea)], 



(22) 



where we have shifted the zero of energy to coincide with 
the bottom of the band, i.e. w(0) = 0. 

We now wish to make an explicit comparison between 
the full response, x{ x ii x 2,t) for this model, and the av- 
eraged response, x(x, i). These quantities can be readily 
computed in one-electron Schrodinger theory, but to keep 
the discussion general we will take a Green's function ap- 
proach. The Fermion field operator for ths system is 



(23) 
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where Ck annihilates an electron with momentum k. 
Since we have only one electron, there is no need to spec- 
ify the spin. In its ground state, the electron is at rest, 
i.e. 



\g> 



4|o> 



(24) 



FIG. 5: x( a; i! x 2, t) averaged over all source locations, 2:2- 
This is the image that would result by inverting IXS data. 
Note the close similarity to Fig 4a. 



Plugging these definitions into the expression for the 
propagator, eq. 1, we find 



5 



X(xi,x 2 ,t) = - [^0(^1)^0(^2)] 

7T 



^ <P(x 2 - n 2 a)(j>(xi - nia) ■ fn 2 - ni (t) 



(25) 



where 



t = 400 ssec 



fn(t) 



7r / a 



dk cos(kna) sinu)(k)t. 



(26) 



In general, this integral should exclude the value at k = 0, 
however we may include it since sinw(0)i = 0. Note that 
fn(t) = /-«(t). 

The propagator, eq. 25, consists of two factors, which 
can be thought of as propagators for the individual elec- 
tron and hole. The first factor 



(27) 



is independent of time and has a strong dependence on 
X2- We plot this factor in Fig. 2, where we have chosen 
(j = a 1 5. We see that this factor is strongly dependent 
on X2, its maximum possible value being reached when 
X2 resides at the center of a well, i.e. the peak of the 
ground state probability density. 

We now evaluate eq. 25 and examine the properties 
of the full propagator. For this purpose we evaluated 
f n (t), out to n = 5 (fifth neighbors), using the value th 
— 1 eV. The results are shown in Fig. 3. The resulting 
x(xi, X2,t) is shown as a contour plot, in Fig. 4, as a 
function of x\ and t, for the specific values X2 — (the 
source in the center of a well) and x 2 = a/2 (source half 
way between two wells) . In accordance with the physical 
meaning of a propagator, at t = the system appears to 
be "struck" by an instantaneous charged perturbation, 
setting up a disturbance in the electron probability den- 
sity. As time progresses the disturbance propagates away 
from the source location, dispersing through the lattice 
in a manner dictated by the dispersion of the band (eq. 
22). 

When X2 = 0, the disturbance is symmetric around 
X\ = and has a distinct temporal periodicity deter- 
mined by the bandwidth. When X2 — a/2, on the other 
hand, the disturbance is asymmetric and no longer has 
a clear temporal periodicity. These two cases are dis- 
tinctly different, and illustrate how local field effects, 
arising from broken translational symmetry, affect the 
response properties of an atomic-like system. 

Most importantly, in Fig. 4c we compare, on a linear 
plot, the two quantities x( x i,0,t) and xi x i, a/2,t) at 
a time t = 1000 as. Notice that the x 2 = case is 
larger than the X2 — a/2 case by a factor of ~ 200. This 
is partly a result of the prefactor eq. 27, and has the 
physical meaning that the system responds more strongly 
if it is "struck" at a region of higher electron probability 
density. 




FIG. 6: x( r i t) f° r a rea l IXS experiment from a single crystal 
of graphite. The six-fold symmetry is a result of the underly- 
ing honeycomb symmetry of the graphene sheets. This image 
can be thought of as a superposition of two sources, one on 
the center of a C atom in each of the crystalline sub-lattices 
(see text). 



We now consider the averaged quantity, x{x, t), which 
was determined by evaluating eq. 17, i.e. by integrat- 
ing images like those shown in Fig. 4a, b for X2 values 
between —a/2 to a/2. The result is shown in Fig. 5. 
Remarkably, it is visually indistinguishable from Fig. 4a. 

The reason these two plots look the same is that 
x( x i, X2,t) is bigger when X2 coincides with the center of 
a well. So the average, eq. 17, is dominated by its con- 
tribution from X2 ~ 0, resulting in a close resemblance to 
Fig. 4a. In other words, the averaged quantity x( x , t) can 
be thought of as the full response xi x i ,X2,t) for x 2 resid- 
ing at the peak in the electron probability density. This 
model is somewhat idealized, but this effect should occur 
quite generally in inhomogeneous systems, so should aid 
interpretation of images generated by IXS inversion. 



V. EXAMPLE: SINGLE CRYSTAL GRAPHITE 

It is enlightening to consider a real example. Graphite 
is a layered structure comprising an abab- type stacking of 
graphene sheets. The sheets have a honeycomb structure 
with two distinct C sites, which are related by a transla- 
tion combined with a 60° rotation of the coordinates. 

Fig. 6 shows a translationally averaged x( r > t), at time 
t = 400as, acquired by inversion of IXS data from a sin- 
gle crystal (not HOPG) of graphite. The images show 
a pronounced six-fold symmetry, which is result of the 
underlying hexagonal symmetry of the lattice. The ex- 
periments were done at Sector 9 XOR of the Advanced 
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Photon Source. The details of this experiment are being 
published in two separate articles [1, For the current 
purposes, we wish only to point out the following. 

The ground state electron density of graphite is maxi- 
mum at the center of a C atom. In accordance with con- 
clusion of Section IV, we expect that Fig. 6 will resemble 
the full x( r i; r 2, t), for the specific case of r2 centered on 
an C atomic site. 

Unfortunately, there are two C atoms in a unit cell. 
They are distinct but related by a simple symmetry op- 
eration, so should have the same electron density. Where 
should one visualize the source residing in this case? Fur- 
ther, the surroundings of each C atom has a three-fold, 
rotational symmetry. Why does Fig. 6 exhibit six-fold 
symmetry? 

The point here is that, by eq. 17, the image Fig. 6 is 
coherent superposition of sources at the two distinct C 
locations. Because they have the same density, they con- 
tribute equally to the average in eq. 17. Further, while 
each C atom locally has only three-fold symmetry, one is 
rotated 60° from the other. The superposition therefore 
exhibits six-fold symmetry, as observed in Fig. 6. Be- 
cause we know that, apart from the 60° rotation, the two 
C atoms are the same, one can in principle reconstruct 
the response of an isolated C atom from Fig. 6. The 
important point for now is that eq. 17 provides a useful 
framework for conceptualizing the meaning of inverted 
IXS data. 

It is worth mentioning that it may be possible to ex- 
tend IXS inversion methods to image more than just an 
average response. Many years ago it was shown jl(| [H| 
that by setting up a standing wave in a material - for 



example by hitting a Bragg reflection - one can mea- 
sure a sub-class of the off-diagonal components of the 
susceptibility. Under standing wave conditions, the ini- 
tial photon will be in a superposition of two different 
momentum eigenstates, in which case the IXS cross sec- 
tion will contain contributions from x(fci, ^2, w), where 
k\ =/= tt2- In principle this method can be used to "an- 
chor the problem" in space, and potentially reconstruct 
the complete x( x ii x ii t)- Whether it is experimentally 
possible to measure a sufficiently large set of these off- 
diagonal terms is a subject of current investigation. 

In summary, we have shown that IXS imaging, as pre- 
scribed in ref. |]|, provides images of the charge propa- 
gator, x(xi, X2,t), that are averaged over all source loca- 
tions, X2- No such average is done in time, however, so 
this method is still a direct probe of electron dynamics 
in the attosecond regime. The size of the propagator is 
much larger when X2 resides in a region of high density, so 
one can think of the averaged images as corresponding to 
the complete response, for x-i residing at the peak of the 
electron density. Standing wave methods, in which the 
incident photon is placed in a superposition of momenta 
may extend the applicability of the technique. 
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